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Abstract 

Zero-Forcing (ZF) has been considered as one of the potential practical precoding and detection method for 
massive MIMO systems. One of the most important advantages of massive MIMO is the capability of supporting a 
large number of users in the same time-frequency resource, which requires much larger dimensions of matrix 
inversion for ZF than conventional multi-user MIMO systems. In this case, Neumann Series (NS) has been 
considered for the Matrix Inversion Approximation (MIA), because of its suitability for massive MIMO systems and 
its advantages in hardware implementation. The performance-complexity trade-off and the hardware implementation 
of NS-based MIA in massive MIMO systems have been discussed. In this paper, we analyze the effects of the 
ratio of the number of massive MIMO antennas to the number of users on the performance of NS-based MIA. 

In addition, we derive the approximation error estimation formulas for different practical numbers of terms of 
NS-based MIA. These results could offer useful guidelines for practical massive MIMO systems. 

I. Introduction 

Massive Multiple-Input Multiple-Output (MIMO) systems were firstly introdueed in [[T]|, and have drawn 
great interest form both aeademia and industry. In sueh systems, eaeh Base Station (BS) is equipped with 
dozens to hundreds of antennas to serve tens of users in the same time-frequeney resouree. Therefore, they 
can achieve significantly higher spatial multiplexing gains than conventional multi-user MIMO systems, 
which offers one of the most important advantages of massive MIMO systems, the potential capability to 
offer linear capacity growth without increasing power or bandwidth IfTl-Pni. 
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It has been shown that, for massive MIMO systems where the number of antennas M, e.g., M = 128, 
is much larger than the number of served users K, e.g., K = 16, [^, [0, Zero-Forcing (ZF) precoding 
and detection can achieve performance very close to the channel capacity for the downlink and uplink 
respectively [|3. As a result, ZF has been considered as one of the potential practical precoding and 
detection method for massive MIMO systems |I2|, Il4l-ll6l. 

For the hardware implementation of ZF, despite of the very large number of M, the main complexity 
is the inverse of a A" x iF matrix 0, [|7]|, [jSl. Unfortunately, for massive MIMO systems, although K is 
much smaller than M, it is still much larger than conventional multi-user MIMO systems. As a result, in 
this case, the computation of the exact inversion of the K x K matrix could result in very high complexity 
[|8l , which may cause large processing delay so that the demands of the channel coherence time is not 
met. Due to this reason, Neumann Series (NS) has been considered to carry out the Matrix Inversion 
Approximation (MIA), because it is well suited for massive MIMO systems and it is advantageous for 
hardware implementation dH, [|7I|, dHlI. 

Despite of the advantages, some potential application issues of the NS-based MIA have also been 
identified. Firstly, for a finite M/K ratio, the NS may not converge, resulting the failure of the algorithm 
0 , [TTl . What M/K ratio could offer high convergence probability is still not clear. Secondly, for the NS- 
based MIA to achieve good performance with quick convergence, the KxK matrix needs to be diagonally 
dominant [[H. In order to satisfy this condition, M K is required [[21, [[8l. Similarly, what M/K ratio 
could provide high probability of diagonally dominant is also not clear. Moreover, with a larger number 
of terms, the NS-based MIA offers closer performance to the exact inversion Q, [I8l|. However, the larger 
number of terms results in more processing cycles. Hence, for practical hardware implementation, the 
number of terms cannot be very large. Although the approximation error analysis was carried out and a 
residual error upper bound of the NS-based MIA was derived IIH, the approximation error analysis with 
high accuracy has not been derived. 

In this paper, we address the three problems listed above. Specifically, we firstly derived a M/K ratio 
condition that offers high convergence probability. Then, we derived another M/K ratio condition that 
provides high probability for the KxK matrix to be diagonally dominant. Finally, we carry out the 
approximation error analysis with high accuracy for practical numbers of terms for the NS-based MIA in 
hardware implementation. 

The remainder of this paper is organized as follows. In Section [III the basis of the NS-based MIA in 
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massive MIMO systems is briefly reviewed. The M/K ratio eondition that provides high eonvergenee 
probability is derived in Seetion Uni Then, another M/K ratio condition that offers high diagonally 
dominant probability for the K x K matrix is derived in Section UVl In Section |Vl the approximation 
error analysis with high accuracy for practical numbers of terms for the NS-based MIA is carried out. 
Finally, after a discussion in Section |Vll conclusions are drawn in Section IVIIi 

II. Basis of NS-Based MIA in Massive MIMO Systems 

Consider a massive MIMO wireless system where the BS is equipped with M antennas to serve K 
single-antenna users in the same time-frequency resource. Then, for the uplink, the M x K channel matrix 
is represented by H = [hmk], where h^k denotes the channel coefficient between the mth antenna and 
the kth user, with m = 1,..., M, and k = 1,..., K. Similarly to [|2l, [171, [|8l, the analysis in this paper 
assumes that the hmk elements are in uncorrelated Rayleigh flat fading, i.e., independent and identically 
distributed (i.i.d.) zero-mean unit-variance complex Gaussian variables. Note that, for the Time-Division 
Duplexing (TDD) mode, due to the channel reciprocity, the downlink has the same channel matrix H as 
the uplink, as long as the transmission duration is within the channel coherence time IfTI-lfbl. 

In order to carry out ZF precoding for the downlink or the ZF detection for the uplink, the pseudo-inverse 
of H needs to be calculated [El, [|ll-[[6l, which is written as 

( 1 ) 

Let G = H^H. Then, in dU), despite of the very large number of M, e.g., 256, in massive MIMO 
systems, the main complexity of the hardware implementation lies in the inversion of the K x K matrix 
G [El, 0, 0. To exploit the large spatial multiplexing gains of massive MIMO systems, although much 
smaller than M, the number of K is much larger than conventional multi-user MIMO systems, e.g., 
K = 16. As a result, the complexity of calculating G“^ may be too high for hardware implementation. 
To address this issue, NS has been considered to carry out the MIA, because it is advantageous in hardware 
implementation and it is suitable for massive MIMO systems 0, 0, 0. Specifically, it can be written 
as 

N-l 

G)^'^5^(I^-©G)"©, (2) 

71=0 

where N denotes the number of terms used in the NS, and S is a K x K diagonal matrix. Note that for 
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(l2l) to work, the requirement below has to be satisfied 


lim {Ik - 0G)'' ^ 0^. (3) 

n—>-oo 


Note that G is a eomplex eentral Wishart matrix beeause the elements of H are i.i.d. eomplex Gaussian 
random variables O. Let a = M/K. As K and M grow, as derived in [ITOll . the largest and the smallest 
eigenvalues of G converge respectively to 


-^max (G) —)■ M 



Amin (G) —)■ M 



(4) 


As a result, if © is chosen as m, which is 


© = — 7 ^^- zIk = -^- Ik, 

M{l + a) M + K ' 

then, 

Amax (©G) —)■ 1 + -— 

1 + va 

Amin (©G) —>■ 1 — -— ~~ 1 =- 

1 + V“ 


(5) 


( 6 ) 


Therefore, the eigenvalues of {Ik — ©G) lie approximately in the range of [—2y/a/{l + a), 2^/a/{l + a)] 
lEl, S' Since 2i/a/(l + a) < 1 when a > 1, the convergence of dH) is satisfied with the choice dS]). 
Moreover, when a is very large, 2i/a/(l + a) —)■ 0, which means that dl]) converges very quickly. Hence, 
a small number of in d2l) can offer close performance to the exact inverse. 

Unfortunately, for finite M and K values, the eigenvalues of the product ©G for a particular channel 
realization can lie outside the range of [— 2y^/(l + a), 2^/a/ (1 + a)] |12|, [|7||, which results in the failure 
of dH). To address this issue, an attenuation factor 6 where 0 < 5 < 1 was introduced in [|3, so dH) changes 
to 


© 


5 

M + K 


Ik- 


(7) 


However, the proper choice of 6 is hard to be determined. On the one hand, if 6 is too large, the non¬ 
convergence issue still exists. On the other hand, if 6 is too small, the convergence speed becomes very 
slow, so the number of N needs to be very large to offer a good MIA, increasing the burden of the 
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hardware implementation. 

Instead of ([7]), another © was applied in |171, which achieves a better MIA [j71|. Specifically, G is 
decomposed as 

G = D + E, (8) 

where D is a diagonal matrix including the diagonal elements of G, and E is a hollow matrix including 
the off-diagonal elements of G. Then, 0 is chosen as 

© = D-^ (9) 

To achieve a good MIA with quick convergence, dH) requires that G is a Diagonally Dominant Matrix 
(DDM) Q, mi, i.e., 

\9ii\>^\9ij\,hj = ( 10 ) 

ijA* 

The performance-complexity trade-off and hardware implementation of the NS-based MIA employing dH) 
have been discussed for the downlink and uplink in [|71 and mi respectively. In both cases, the NS-based 
MIA employing dH) was considered as a promising and practical method for massive MIMO systems. As 
a result, the analysis carried out in this paper is based on the choice of dH). 

As mentioned in Section HI there still some issues on the application of dH) for finite M and K values. 
Firstly, it is unclear that what a can offer high convergence probability. Secondly, it is unclear that what a 
can achieve high probability for G to be diagonal dominant. Moreover, more accurate approximation error 
analysis for practical N values is needed. In the next sections, the aforementioned issues are addressed. 

III. Convergence and a 

According to the theory of matrix power series @, for a iT x iT matrix B, the product converges 
to Qk only when the spectral radius of B, denoted by p(B), i.e., the maximum modulus of eigenvalues 
of B, is less than 1. Then, for the choice of dH), a good MIA of dU requires 

p(li,-D-iG) < 1. (II) 

Since the elements of H are i.i.d. zero-mean unit-variance complex Gaussian random variables, when the 
number of M is large, the diagonal elements of D approach to M'Ej{\hkk\^} = M by the law of large 
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numbers |[I1, [|3, flU]. Therefore, the diagonal matrix D ean be replaeed by MIk-, Then, the eondition 
(fTTI) ehanges to 

|M-A(G)| < M ^ 0 < A(G) < 2M. (12) 


As G is a positive-definite matrix IIH, its eigenvalues are all larger than 0 [|9]|. As a result, (fT^ is equivalent 
to 

A^^(G) < 2M. (13) 


Note that G = is a eomplex eentral Wishart matrix |0, and the distribution of Amax(G) is 

provided in IfTTl as 


P (Amax(G) < x) 


CTk {K) km 
CTk (M + K) 


X iFi (M; M + K] -xl) , 


(14) 


where a: is a non-negative number. The eomplex multivariate gamma funetion CTp{a) is defined as 

p 

CTp (a) = 7rP(P-i)/2 JJ r [a - ? + 1], (15) 

i=l 

where p is a positive integer, a is a eomplex-valued number, and r[a] is the gamma funetion. The 
hypergeometrie funetion iFi(M; M + K] —xl) is 


,F,{M-M + K--x1) = J2Y1 

k=0 K 


[Ml C.(-rf) 
[M + K\^ k\ 


( 16 ) 


The details of [M]^ and Ci^{—xT) in (fT^ ean be found in IfTTl . Based on (fT4l) . the probability of (fT3l) ean 
be direetly derived. However, (fT4l) ineludes the summation of infinite terms in (fTb]) whieh has extreme 
eomplexity, so it eannot provide a elosed-form eonvergenee eondition of (fT3l) in terms of a. 

Fortunately, based on (fH), the eondition of (fT3l) ehanges to 


M( 1 + — 

'a 


< 2M. 


(17) 


Based on (fTVl) . a high probability eonvergenee eondition in terms of a is derived as 

a >-i, ss 5.83. (18) 


With (fTSl) . the maximum possible number of K ean be found for a speeifie number of M to aehieve a 
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Fig. 1. The maximum K values for different M values that satisfy GUl 


TABLE I 

Typical M values with their associated maximum values of K and convergence probability values 


M 

64 

128 

256 

512 

K 

10 

21 

43 

87 

Probability of dH) 

0.999 

0.998 

0.995 

0.991 


very high probability of convergence for ([3]). 

Fig. m illustrates the maximum values of K corresponding to M values that vary from 64 to 512 based 
on the convergence condition (flSl) . With these K values, the simulated convergence probability values 
of dU) based on the accurate condition (fTTI) and the approximated condition (fTSl) are shown in Fig. [21 
The results indicate that they provide close probability with (ITTI) being slightly better in massive MIMO 
systems with large M. The results verify that (fTSl) is an acceptable approximation of ([TTI) . Furthermore, 
the results show that the condition (IT^ in terms of a can offer high convergence probability for (|3]). Table 
jl] summarizes the typical M values of massive MIMO systems with their corresponding maximum values 
of K and the convergence probability values of ([3]). Note that ([T^ does not ensure fast convergence of 
(|3]), so a more strict a condition for G being a DDM is studied in the next section. 







































Fig. 2. Convergence probability values of (IHl for the K values in Fig. [T] 


IV. Diagonally Dominant and a 


Let hfc denote the fcth eolumn veetor of the M x K ehannel matrix H. Then, represents the M- 
dimensional channel vector for the kth user. Hence, the elements of the K x K matrix G = is 

calculated as 


9a 11^*11 2 ’ ^ • • •) 

gij = hfhj, j = j 


(19) 


As mentioned in Section Hill the diagonal elements gu approach to M when the number of M is large. 
As a result, the requirement (ITOl) in Section HI] for G being a DDM can be approximated as 


Ai = ^ IAjI < 1, Vi, 

where rij is the normalized correlation coefficient between hj and defined as 




Tij = 


b^b 

iij iij 


0112 


( 20 ) 




M 


( 21 ) 







































Let X = \rij\. Note that the Probability Density Function (PDF) of x was derived in ifT^ as 


f (x) = 2 {M — 1) a: (l — ^ , 0 < x < 1. 


(22) 


Hence, the mean of x is 


E (x) = / x/(x) dx = (M — 1) B (1.5, M — 1), 


(23) 


where B(a, b) with a and b being complex-valued numbers is the beta function defined as 


B(a,6)= / >0. 


(24) 


Although (1231) provides the values of E(x), since the number of K is not large enough, Aj in (l20l) can 
be larger than {K — l)E(x). However, Aj has a high probability being smaller than (K — l)[E(x) + ^(x)] 
where 5{x) denotes the standard deviation of x, which is 


5 (x) = 


(25) 


with 


E(x 2)= / xV(2:)dx = (M-1)B(2,M-1). 


(26) 


Therefore, the condition (l20l) can be approximated as 


{K - 1) [E (x) + 6 (x)] < 1. 


(27) 


Based on (l27l) . a high probability condition for the G matrix being a DDM in terms of a is derived as 


M [E (x) + 5 (x)] 
E (x) + 5 (x) + 1 ’ 


(28) 


With (|2^ . the maximum possible number of K can be found for a specific number of M to achieve a 
very high probability for G being a DDM. 

Fig. [3] shows the maximum values of K corresponding to M values that vary from 64 to 512 based on 
the diagonally dominant condition (l28l) . With these K values, the simulated DDM probability based on 
the definition (fTOl) and the approximated condition (l20l) are illustrated in Fig. IH The results show that they 
achieve close probability in massive MIMO systems with large M. The results verify that (l20l) is a good 
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Fig. 3. The maximum K values for different M values that satisfy l l28b 


TABLE II 

Typical M values with their associated maximum values of K and diagonally dominant probability values 


M 

64 

128 

256 

512 

K 

6 

9 

12 

17 

Probability of (fl^ 

0.990 

0.977 

0.998 

0.999 


approximation of (fT^ . especially when M is very large. Moreover, the results show that the condition 
(1281) in terms of a can offer high DDM probability. Table HI] summarizes the typical M values of massive 
MIMO systems with their corresponding maximum values of K and the diagonally dominant probability 
values of (fTOl) . Note that the DDM condition (|2^ is sufficient for the convergence condition (flSl) and 
leads to quicker convergence, so it is more useful in practice. 

V. Error Analysis 

Based on (jS]) and dH), the NS-based MIA of (O changes to 

Af-l 


n 


(29) 
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Fig. 4. Diagonally dominant probability values for the K values in Fig. [3] 


Note that if the eonvergenee eondition ([3]) is satisfied, is the exaet matrix inverse of G. However, in 
praetiee, the number of N eannot be very large. Otherwise, it would eause exeessive burden for hardware 
implementation. In this ease, residual error resulted from the NS-based MIA G^^ exists. Let the K- 
dimensional veetor s denote the transmitted symbols for the uplink or the downlink. Without loss of 
generality, Eds^P) = 1 is assumed, with k = 1,.. .,K. Let Z = D“^E. Then, the Mean Square Error 
(MSE) of the NS-based MIA G^^ for the uplink is derived as 


€S = e{||(G-‘-G;;‘) H“Hs||d 

oo 

Z^^(-Z)"D-^H^Hs 


= E 

n=0 

= e||1z^g-^Gs||^ 


= E<^ llz^s 


= E|Tr Z^ss^(Z^)” I 
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= E<^Tr 


H(2iV)Hz7vj| 

= Tr|E [ss^] E (Z^)“Z^ | 
= Tr|li^E (Z^)“Z^ I 

(z^)“z^ 


= E<^Tr 



Note that for the downlink ease, the MSE result is 


(30) 


4 = e{||sT(G„‘-G„‘)H>‘H||^} 

= e{||s’'z'"||^} 

= e{||Z'»||p}, (31) 


whieh is the same as (1^ . Henee, eAi is used instead of and in this seetion below. Sinee 
ean be interpreted as the power of the residual interferenee of ZF preeoding or deteetion, the average 
Signal-to-Interferenee Ratio (SIR) for eaeh user is ealeulated as 


In = 


NP 

K 

ejv 

K 


K 


In lISl, the MSE eAr in (1^ and (ISTl) is upper bounded as 


(32) 


€-N < 


{K^-K), 


2M(M + 1) 


(M - 1) (M - 2) (M - 3) (M - 4) 


N 


(33) 


with M > A. Unfortunately, (l3^ is a very loose upper bound, resulting in a very loose lower bound of 
7 Ar in (l3^ . In Fig. |5l the exaet SIR values and the lower bound values are eompared with M = 128 for 
different K and M values. The results show substantial differenees when > 1, whieh eannot provide 
suffieient insight for the residual error of the NS-based MIA for > 1. Due to this reason, we seek to 
derive a more aeeurate approximation of eTv in this seetion below. 

When M is large, beeause D ean be approximated as Ml^ as mentioned in Seetion |nll aeeording to 
(fT9l) and (l2T]) . the elements of Z is approximated as 


0 


i = 1, 


K. 


= Zji 


h^h, 

1 J 

~ M ~ 


Uj, j = j 


(34) 
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Fig. 5. Comparison between the exact SIR values and the lower bound values resulting from 03l l with M = 128 for different K and N 
values. 


As a result, the PDF of x = \zij\ ean be approximated as (l22l) . Then, a more aeeurate approximation of 
7 Ar can be derived based on (l22l) . 

When = 1, the MSB cat in (1^ and (l3TI) changes to 

K K 

e, = \\Z\\l = J2 \^ijf-K{K-l)E{x^). (35) 

i=l j=l,jjti 

Since E(x^) has been derived as (|2^ . the term ei in (1351) is rewritten as 

ei^ K{K -1) B2,m, (36) 


where Ba,M is defined as 


B.,m = (M-l)B(a,M-l). 


When = 2, the MSB cat in (l30l) and (l3TI) changes to 


(37) 


£2 = 


(38) 
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where the elements in Y = is 



k=l,k^i 


(39) 


K 


Vij = E ^ikZjk, j = j 


k=l,k^ij 


Note that ||Z^||| ean be written as a summation of polynomial terms, whieh ean be elassified into three 
eategories. The first category includes {K — 1)K terms of with i ^ k. The second category includes 
{K — 2){K — 1)K terms of with i 7 ^ fc 7 ^ /, as well as {K — 2){K— 1)K terms of l-SjfcPkjfcP 

with i ^ j ^ k. Hence, the total number of terms for the second category is 2{K — 2){K — 1)K. Finally, 
the third category includes {K — 3){K — 2){K — 1)K terms of ZikZ*^z*iZji with i ^ j ^ k ^ 1. Because 
the elements of H are i.i.d zero-mean unit-variance complex Gaussian random variables, based on (l34l) . 
the elements of Zij are i.i.d. zero-mean random variables. As a result, the terms of the third category 
are also i.i.d. zero-mean random variable. Therefore, the sum of the terms of the third category can be 
approximated as zero. For the terms of the first category, the mean can be calculated based on (l22l) as 



(40) 


Similarly, the mean of the terms of the second category is can be approximated as 



(41) 


Due to (l40l) and (l4T]) . the term 62 in (l3^ is approximated as 


e2^K{K-l) + 2{K-2){K-1) KBI^. 


(42) 


When N > 2, the MSB cat in (l30l) and (|3TI) can be derived with the similar method applied by Y = 2. 
The results of Y = 3 and Y = 4 are directly provided below as 


63 ^ (Y - 2) (Y - 1) Y (5Y - 8) Yj^ 

+ (2Y — 3) (Y — 1) KB3^mB2,m, 


(43) 


and 


64^(2Y-3) {K-l)KBl^ 
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Fig. 6. Comparison between the exact and estimated SIR values with M = 128 for different K and N values. 

+ {K-2){K-lfK‘^Bl^. (44) 

With the estimated residual error formulas (l3^ . (I42l) - (l44l) . the estimated SIR formulas ean be easily 
derived aeeording to (l32l) . Fig. [MH] eompare the exaet and estimated SIR values for different K and N 
values, with M = 128, M = 256, and M = 512 respeetively. The results show that the estimated SIR 
values are very elose to the exaet SIR values, whieh verifies the high aeeuraey of SIR estimation formulas 
based on (1^ . (I42l)-(l44l). 


VI. Discussions 

In massive MIMO systems, a is eommonly eonsidered to be very large to offer good performanee [|2l, 
[14| , e.g., a > 10. Henee, the eonvergenee eondition (fT^ . i.e., a > 5.83, derived in Seetion|ni]is generally 
satisfied for massive MIMO systems. Note that the eonvergenee probability values provided in Fig. [T] and 
Table HI whieh are already elose to 1, eorrespond to the smallest a values that satisfy (fT^ . Henee, the 
eonvergenee probability values for massive MIMO systems are not lower than the values provided in Fig. 
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Fig. 7. Comparison between the exact and estimated SIR values with M = 256 for different K and N values. 


[Hand Table HI Therefore, the eonvergenee of NS-based MIA is guaranteed so that it is a valid method for 
massive MIMO systems, and its accuracy can be improved by increasing N. 

As mentioned at the end of Section [Till the convergence condition ([T^ does not guarantee quick 
convergence of ([3]). With the diagonally dominant condition (l28l) derived in Section jlVl however, the NS- 
based MIA can achieve good accuracy with quick convergence, i.e., a small N can offer a sufficiently good 
MIA. Otherwise, with the same N value, violating (1281) results in performance loss for the ZF decoding 
or detection employing the NS-based MIA. Take the simulation results provided in [O as examples, with 
M = 128 and N = 3, the choice of K = A satisfying (1281) achieves close performance to the exact 
inverse, while the choice of iT = 12 violating (1281) suffers huge performance loss. However, (1281) requires 
very small a values, and a becomes smaller as M increases, which can be seen from Table HVl The strict 
requirement of a may reduce the spatial multiplexing advantage of massive MIMO systems, i.e., at most 
K = 17 users can be served by M = 512 antennas. To relieve this issue, one comprised choice is to apply 
an a slightly higher than (|2^ with slightly larger N of the NS-based MIA, depending on the hardware 
capability. 
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Fig. 8. Comparison between the exact and estimated SIR values with M = 512 for different K and N values. 


The SIR discussed in Section |V] reflects the performance error floor for ZF precoding or detection 
employing practical NS-based MIA in massive MIMO systems. The performance error floor decides the 
best performance that the ZF precoding or detection employing the NS-based MIA can achieves. As 
a result, with M, K, and N, the best achievable performance can be easily estimated based on (1^ . 
(I42l) - (l4^ . In addition, since larger N causes higher hardware implementation complexity, with M, K, 
and the target performance, the smallest choice of N that can offer sufficiently good performance can be 
determined to relieve the complexity. Note that a revised form of © was provided in [|3 as 


Af-l 


5^(1^-©or© 


n=0 

L-1 


H [Ik + (Ik - 0 G)'] ©, 
1=0 


(45) 


where L is a positive integer with N = 2^. Hence, L = 1, L = 2, and L = 3 of the alternative expression 
(1451) correspond to = 2, = 4, and A^ = 8 of the regular expression © respectively. As a result, with 

the alternative expression (145]) . after the choice of A^ = 4, the NS-based MIA with the choice of A^ = 8 
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can be quickly calculated. Therefore, if the choice of iV = 4 is not good enough based on the estimation 
formula dH]), the choice of iV = 8 can be directly selected based on (1431) . Furthermore, note that the 
complexity of the NS-based MIA with the choice of > 3 is considered to be O(iT^) in |[8l, which 
loses the complexity advantage over the exact matrix inverse of 0{K^). In fact, however, the NS-based 
MIA can be implemented as a series of cascaded matched filter so that the complexity can be reduced to 
0(A'^), as discussed in In this way, the NS-based MIA still has the complexity advantage over the 
exact inverse even with the choice of = 8. 


VII. Conclusions 

In this paper, three issues related to the practical application of the NS-based MIA in massive MIMO 
systems are addressed. Firstly, a > 5.83 as in (fT^ is offered for the NS-based MIA to achieve very high 
convergence probability. In other words, with the number of BS antennas M, the maximum number of 
served users K for the NS-based MIA to be a valid method in massive MIMO systems can be determined. 
Then, a tighter condition (1281) is provided for G to be a DDM in very high probability, resulting in a 
good NS-based MIA with a small number of N. This means that given the number of BS antennas M, 
the maximum number of served users K for the NS-based MIA to achieve good performance and quick 
convergence for ZF decoding or detection cab be determined. Finally, by approximation error analysis, 
residual error estimation formulas (1^ . (I42l)-(l4^ with very high accuracy are derived for practical N 
values, which can be applied to estimate the error floor caused by the NS-based MIA. Thus, given the 
number of BS antennas M, the number of served users K, and the number of terms employed by the 
NS-based MIA N, highly accurate estimation of the SIR caused by the NS-based MIA can be obtained. 
These results offer useful guidelines for practical application of the NS-based MIA in massive MIMO 
systems. 


References 

[1] T. L. Marzetta, “Noncooperative Cellular Wireless with Unlimited Numbers of Base Station Antennas,” IEEE Trans. Wireless Commun., 
vol. 9, no. 11, pp. 3590-3600, Nov. 2010. 

[2] R Rusek, D. Persson, B. K. Lau, E. G. Larsson, T. L. Marzetta, O. Edfors, and F. Tufvesson, “Scaling up MIMO: Opportunities and 
Challenges with Very Large Arrays,” IEEE Signal Process. Mag., vol. 30, no. 1, pp. 40-46, Jan. 2013. 

[3] E. G. Larsson, F. Tufvesson, O. Edfors, and T. L. Marzetta, “Massive MIMO for Next Generation Wireless Systems,” IEEE Commun. 
Mag., vol. 52, no. 2, pp. 186-195, Feb. 2014. 

[4] J. Hoydis, S. Brink, and M. Debbah, “Massive MIMO in the UL/DL of Cellular Networks: How Many Antennas Do We Need?” IEEE 
Sel. Areas Commun., vol. 31, no. 2, pp. 160-171, Feb. 2013. 


19 


[5] C. Shepard, H. Yu, N. Anand, L. E. Li, T. Marzetta, R. Yang, and L. Zhong, “Argos: Practical Many-Antenna Base Stations,” in Proc. 
MobiCom 12, Istanbul, Turkey, Aug. 2012. 

[6] H. Yang and T. L. Marzetta, “Performance of Conjugate and Zero-Forcing Beamforming in Large-scale Antenna Systems,” IEEE J. 
Sel. Areas Commun., vol. 31, no. 2, pp. 172-179, Feb. 2013. 

[7] H. Prabhu, 1. Rodrigues, O. Edfors, and F. Rusek, “Approximative Matrix Inverse Computations for Very-large MIMO and Applications 
to Linear Pre-coding Systems,” in Proc. IEEE WCNC 13, Shanghai, China, Apr. 2013. 

[8] M. Wu, B. Yin, G. Wang, C. Dick, 1. Cavallaro, and C. Studer, “Large-Scale MIMO Detection for 3GPP LTE: Algorithms and FPGA 
Implementations,” IEEE J. Sel. Topics Signal Process., vol. 8, no. 5, pp. 916-929, Oct. 2014. 

[9] R. A. Horn and C. R. Johnson, Matrix Analysis. Cambridge Univeristy Press, 1990. 

[10] A. Edelman, “Eigenvalues and Condition Numbers of Random Matrices,” Ph.D. dissertation, MIT, 1989. 

[11] T. Ratnarajah, R. Vaillancourt, and M. Alvo, “Eigenvalues and Condition Numbers of Complex Random Matrices,” SIAM Journal on 
Matrix Analysis and Applications, vol. 26, no. 2, pp. 441-456, 2004. 

[12] D. Zhu, B. Li, and P. Liang. (2014) Normalized Volume of Hyperball in Complex Grassmann Manifold and Its Application in 
Large-Scale MU-MIMO Communication Systems. arXiv: 1402.4543. [Online]. Available: http://arxiv.org 


